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1.  Introduction 


The  effect  of  strain  rate  on  ductile  fracture  is  one  of  the  least  understood  phenomena  in 
modem  fracture  mechanics.  At  the  same  time,  information  on  ductile  fracture  under  dynamic 
loading  is  very  important  for  reliable  prediction  of  fracture  and  fragmentation  of  high- 
consequence  structures.  High-consequence  structures  are  understood  here  as  components  of 
turbofan  engines  for  fixed-wing  aircraft  and/or  rotorcraft  dynamic  components.  The  high 
strain  rate  experiments  to  fracture  are  conducted  with  equipment  available  at  the  Institute  of 
Fundamental  Technical  Research. 

Mechanical  and  thermomechanical  behavior  of  materials  (DH-36  steel)  applied  in 
airplanes  and  helicopters,  subjected  to  compression,  tensile  and  simple  shear  tests  in  wide 
spectrum  of  the  strain  rates  (10~3  s'1  4-  104  s'1)  are  investigated.  The  tensile,  cycling  and  shear 
properties  seams  to  be  very  important,  since  the  aircraft  materials  are  subjected  to  such  a 
loading  during  their  application.  The  simple  shear  test  of  the  DH-36  steel  sheets  in  the  case  of 
quasi-static  case  as  well  as  in  the  high  strain  rates  is  performed.  In  this  kind  of  test  (proposed 
by  W.K.  Nowacki  [1])  a  considerable  homogeneity  of  the  permanent  strain  field  at  finite 
deformations  over  the  total  length  of  the  specimens  is  observed  in  experiments  of  simple 
shear  and  in  the  results  of  simulation  -  cf.  [2],  The  proposed  method  is  the  only  known  test 
providing,  in  the  case  of  a  thin  sheet,  homogeneous  stress  and  strain  fields  in  both  the 
dynamic  and  static  tests.  Simple  shear  test  is  particularly  attractive,  since  the  application  of 
this  type  of  loading  path  can  result  in  large  strains  without  the  occurrence  of  plastic 
instability.  The  shear  tests  on  small,  flat  specimens  using  are  described  in  the  Chapter  2. The 
mechanical  characteristics  and  the  temperature  distributions  of  the  sample  surface  obtained  by 
a  thermovision  camera  during  the  deformation  processes  are  recorded,  processed  and  stored  in 
the  computerized  registration  system.  The  computer  results  are  available  for  further  analysis; 
its  later  reproduction  in  the  most  convenient  is  being  possible.  Investigation  of  the  initiation 
and  further  development  of  fracture  locus  at  various  strain  rates  using  the  infrared  technique 
are  performed.  Identification  of  the  beginning  of  plastic  deformation,  based  upon  the 
temperature  measurements,  in  the  materials  subjected  to  subsequent  shear  deformation  are 
carried  out. 
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The  objective  of  the  Chapter  3  is  to  discuss  very  efficient  procedure  of  the  numerical 
investigation  of  localized  fracture  in  inelastic  solids  generated  by  impact-loaded  adiabatic 
processes.  Particular  attention  is  focused  on  the  proper  description  of  a  ductile  mode  of 
fracture  propagating  along  the  shear  band  for  high  impact  velocities.  This  procedure  of 
investigation  is  based  on  utilization  the  finite  element  method  and  ABAQUS  system  for 
regularized  thermo-elasto-viscoplastic  constitutive  model  of  damaged  material.  A  general 
constitutive  model  of  thermo-elasto-viscoplastic  damaged  polycrystalline  solids  with  a  finite 
set  of  internal  state  variables  is  used.  The  set  of  internal  state  variables  consists  of  two  scalars, 
namely  equivalent  inelastic  deformation  and  volume  fraction  porosity.  The  equivalent 
inelastic  deformation  can  describe  the  dissipation  effects  generated  by  viscoplastic  flow 
phenomena  and  the  volume  fraction  porosity  takes  into  account  the  microdamage  evolution 
effects.  The  relaxation  time  is  used  as  a  regularization  parameter.  Fracture  criterion  based  on 
the  evolution  of  microdamage  is  assumed. 

As  a  numerical  example  we  consider  dynamic  simple  shearing  and  localized  fracture 
in  thin  plate.  The  impact  loading  is  simulated  by  a  velocity  boundary  condition  which  are  the 
results  of  dynamic  contact  problem.  The  separation  of  the  projectile  from  the  specimen, 
resulting  from  wave  reflections  within  the  projectile  and  the  specimen,  occurs  in  the 
phenomenon. 

A  thin  shear  band  region  of  finite  width  which  undergoes  significant  deformation  and 
temperature  rise  has  been  determined.  Its  evolution  until  occurrence  of  final  fracture  has  been 
simulated.  Shear  band  advance,  shear  band  velocity  and  the  development  of  the  temperature 
field  as  a  function  of  time  have  been  determined.  Qualitative  comparison  of  numerical  results 
with  experimental  observation  data  has  been  presented.  The  numerical  results  obtained  have 
proven  the  usefulness  of  the  thermo-elasto-viscoplastic  theory  in  the  investigation  of  dynamic 
shear  band  propagations  and  localized  fracture. 

In  Chapter  3  an  application  of  certain  particular  model  of  viscoplasticity,  in  which  a 
multiscale  mechanism  of  shear  banding  is  accounted  for,  is  presented.  Also  a  possible  method 
of  identification  of  this  model  is  discussed.  The  identification  of  the  proposed  model  can  be 
based  on  the  experimental  results  presented  in  [15]. The  presented  model  can  be  applied  for 
the  analysis  of  strain  rate  effect  on  ductile  fracture  of  DH-36  steel  with  use  of  the  original 
methodology  developed  by  T.  Wierzbicki  and  coworkers,  cf.  eg.  [16]. 
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2.  New  shear  fracture  tests  on  small,  flat  specimens  using  a  unique  apparatus  developed 
at  IPPT 


In  the  static  investigations  of  thin  sheets  we  can  distinguish  two  approaches,  based  on 
different  principles  of  the  shear  device.  The  first  type  of  shear  device  of  the  one-shear  zone 
only,  proposed  by  G’Sell  et  al.  [3],  has  been  performed  for  the  investigation  of  thin  polymer 
sheets  in  the  conditions  of  static  shear.  The  principle  of  this  construction  has  been  adopted  for 
testing  of  the  steel  sheets,  Tourabi  et  al.  [4]  and  Rauch  and  G’Sell  [5],  cyclic  shear  of  metallic 
sheets,  Wack  et  Tourabi  [6]  and  for  the  shape-memory  alloy,  Manach  and  Favier  [7].  In  the 
second  type  of  device  for  the  static  tests  (see  Fig.  1),  Yoskida  &  Myauchi  [8],  the  Rouble- 
shear”  specimen  is  applied.  It  is  a  specimen  with  two  zones  of  shear.  The  fundamental 
attribute  of  this  device  is  its  symmetry  with  respect  to  the  direction  of  loading.  Both  the  above 
experimental  techniques  have  their  advantages  and  the  inconveniences.  In  the  first  method, 
proposed  by  G’Sell,  one  can  obtain  a  uniform  zone  of  the  deformation  field.  First  of  all,  it  is 
possible  in  this  method  to  select  one  direction  in  which  we  can  examine  the  behavior  of 
material.  The  symmetry  of  loading  may  need  special  very  sophisticated  carriages.  The  second 
method,  proposed  by  Myauchi,  is  the  simplest  in  use.  In  this  method,  the  rigorous  conditions 
of  symmetry  may  be  fulfilled.  The  exact  fixation  in  grips  must  be  assured.  The  inconvenience 
of  this  method  is  the  fact  that  the  homogeneity  of  the  deformation  field  is  not  very  good, 
compared  with  the  first  method.  The  free  boundaries  of  the  specimen  have  the  influence  for 
the  field  of  deformations  and  stresses. 

In  the  case  of  dynamic  loading,  using  the  Kolsky  bar,  also  called  a  Split  Hopkinson  Pressure 
Bar  (SHPB),  the  scheme  of  the  shear  device  (see  Fig.  1)  proposed  by  Myauchi,  is  more 
convenient  for  this  kind  of  test.  Regarding  the  miniature  of  the  shear  device,  the  scheme  with 
horizontal  and  vertical  carriages,  proposed  by  G’Sell,  is  practically  unfeasible.  A  new  shear 
device  was  used  by  Gary  and  Nowacki  [1]  to  perform  tests  under  high  strain  rates  on 
specimens  having  the  form  of  slabs,  such  as  metal  sheets.  The  loading  and  the  displacements 
of  this  device  are  controlled  by  a  Split  Hopkinson  Pressure  Bar  acting  in  compression.  The 
role  of  the  special  device  is  to  transform  the  compression  into  a  plane  shear  (as  in  the 
Myauchi  system).  The  dynamic  simple  shear  tests  so  obtained  is  the  only  known  method 
allowing  to  obtain  a  very  good  homogeneity  of  the  residual  strain  field  over  the  total  length  of 
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the  specimen,  without  localization  of  the  deformations,  as  in  the  case  of  torsion  of  thin-walled 
tubes,  Marchand  and  Duffy  [9]. 


Figure  1.  Scheme  of  shear  testing;  two  shear  zones 

The  shear  device  for  investigations  at  high  strain  rates  consists  of  two  coaxial  cylindrical  parts, 
namely  the  external  tubular  part  and  the  internal  massive  part  (see  Fig.  2).  Both  cylinders  are 
divided  into  two  symmetrical  parts,  and  between  them  the  sheet  to  be  tested  is  fixed. 


Figure  2.  Shear  device 


Two  boundaries  of  the  specimen  between  the  internal  and  external  parts  of  the  device  are  in 
plane  shear  when  these  cylinders  move  axially  one  toward  the  other.  Each  shear  zone  before 
test  is  rectangular  and  becomes  very  close  to  a  parallelogram  having  a  constant  length  and  a 
constant  height  (see  Fig.  1  and  Fig.  3). 
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time  t  =  0  ^ 


incident  bar 


time  t,  >  0 


transmitter  bar 


time  t2  >  time  t1  "°  a0  and  h0 - >  constant 


Figure  3.  Scheme  of  simple  shear.  Loading  and  displacement  are 
controlled  by  a  SHPB  acting  in  compression 


The  greatest  double-shear  specimen  height  allowed  by  the  shear  device  is  l0  =  30  mm.  The 
width  value  a„  is  chosen  with  a  view  to  satisfying  two  requirements:  first  the  ratio  ajd  of  the 
width  to  the  thickness  has  to  be  small  enough  to  avoid  buckling  of  the  sheared  zone  (a„/d  < 
10)  and  secondly,  the  ratio  a(/l„  of  the  width  to  the  height  has  to  be  sufficiently  small  to 
minimize  the  error  due  to  the  no  homogeneity  of  the  shear  stress  and  strain  at  the  two  ends  of 
the  sample  (a0/l„<  l/\  0),  Tourabi  et  al.  [4]  and  Gary  and  Nowacki  [1],  In  our  case,  the  material 
sheet  has  a  nominal  thickness  of  0.64  mm  (usually  specimens  can  have  different  thickness).  In 
our  experimental  investigations  we  used  specimens  made  of  DH-36  steel.  We  apply  in  our 
quasi-static  tests  the  ratio  aflg  =  1/10  (the  width  of  the  shear  zone  a0  is  constant  and  equal  to 
3mm;  ajd  =  4.7).  In  the  case  of  dynamic  test  at  high  strain  rates  we  use  specimens  of  the 
strain  gauge  length  20mm  [10].  Design  of  the  specimen  for  dynamic  testing  in  simple  shear  is 
presented  on  the  Fig.  4. 

It  is  important  to  mention  that  we  use  our  shear  device  in  a  large  range  of  strain  rates:  in  the 
case  of  quasi-static  tests  from  10'  s' ,  to  10  s'  in  the  case  of  dynamic  tests  at  high  strain 
rates.  Using  the  same  shear  device  in  all  experiments,  we  can  obtain  a  conviction  of  the 
reliability  of  our  experimental  results. 


7 


Figure  4.  Specimen  used  during  double  shear  tests.  On  the  figure  b)  specimen  after  and 

before  simple  shear  deformation 

2.1.  Quasi-static  tests  in  simple  shear.  Deformation  processes  always  modify  the 
temperature  fields  in  materials  while  the  temperature  variations  influence  the  stress-strain 
characteristic.  Strong  effects  of  the  thermomechanical  coupling  are  observed  during  the 
process  of  the  simple  shear  test.  During  this  test  the  first  invariant  of  stress  tensor  Ao\\  =  0, 
thus,  in  elastic  range  of  deformation  temperature  does  not  change  [11],  The  following  plastic 
deformation  is  accompanied  by  increase  in  temperature.  It  is  caused  by  transformation  of  a 
significant  part  of  the  plastic  energy  into  heat  [12], 

The  goal  of  our  investigation  was  obtaining  the  mechanical  curves  as  well  as  the 
temperature  distributions  in  the  shear  areas  and  comparison  them  in  order  to  analyse  the 
processes  of  materials  deformation.  The  samples  were  placed  in  a  specially  designed  grip  for 
quasi-static  tests,  transforming  the  compression  into  a  simple  shear  process  (see  Fig.  5).  The 
grip  was  fitted  in  the  testing  machine,  working  in  compression  mode.  A  change  of 
temperature  of  the  surface  of  two  shear  paths  has  been  observed  during  the  test  using  an 
infrared  camera. 
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shear  zones 


Figure  5.  Photograph  of  shear  test  device 

The  tests  were  performed  at  room  temperature  (295  K),  with  the  rates  of  deformation 
equal  to  10 3  s'1,  10'1  s'1  and  10°  s'1. 

During  deformation,  the  load  vs.  crosshead  displacement  and  the  distribution  of  infrared 
radiation,  emitted  by  the  specimen,  were  continuously  registered.  The  distribution  was 
measured  using  the  thermovision  camera  coupled  with  a  computer  system  of  data  acquisition 
and  conversion. 

The  mechanical  characteristics  made  it  possible  to  derive  the  stress-strain  relations.  The 
distribution  of  intensity  of  the  infrared  radiation  registered  in  the  same  time  allowed  us  to 
reconstruct  thermal  pictures  (Fig.  6)  of  the  samples  and  in  this  way  to  find  the  temperature 
changes  of  the  sample  subjected  to  shear  test.  The  increments  of  the  sample’s  temperature 
were  calculated  as  an  average  value  from  the  central  segments  along  the  sheared  zones.  This 
approach  let  us  to  compare  the  mechanical  and  the  thermal  characteristics,  and  then  to  find 
temperature-stress  relations  AT{pu). 
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Figure  6.  Example  of  thermogram  obtained  during  simple  shear  test  of  DH-36  steel 

The  mean-square  error  of  temperature  evaluation  was  about  0.2  4-  0.5  K.  In  order  to  ensure 
higher  and  more  homogeneous  emissivity,  the  surfaces  of  the  shear  zones  of  the  steel  sheet 
specimens  were  blackened  with  a  carbon  powder. 

As  results  of  the  study  we  have  obtained  stress-strain  and  temperature-strain  characteristics. 
Example  of  these  stress-strain  relations  obtained  for  DH-36  steel  deformed  with  shear  rate  equal 
to  10~3  s'1,  10'1  s'1  and  10°  s"1  are  presented  in  Figs.  7a,  b  and  c  respectively  and  temperature- 
strain  relations  are  given  in  Figs.  8a  and  b,  for  strain  rates  10'1  s'1  and  10°  s'1. 
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Figure  7a.  Stress-strain  relation  obtained  during  simple  shear  of  DH-36  steel  with  shear 

3  1 

rate  equal  10'  s' ;  2,3,4  -  tests  numbers 
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Shear  stress  (MPa) 


Figure  7b.  Stress-strain  relation  obtained  during  simple  shear  ofDH-36  steel  with  shear 

rate  equal  I  O' !  s'1;  5,7,8  -  tests  numbers 


Figure  7c.  Stress-strain  relation  obtained  during  simple  shear  ofDH-36  steel  with  shear 

rate  equal  10°  s'1;  10,11,12  -  tests  numbers 
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Figure  8a.  Stress-strain  and  mean  temperature  changes  curves  ofDH-36  steel  at  the  rate 

of  deformation  equal  l  O'1  s'1 
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Figure  8b.  Stress-strain  and  mean  temperature  changes  curves  ofDH-36  steel  at  the  rate  of 

deformation  equal  10°  s'1 
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Figure  9  a).  Stress-strain  relation  obtained  during  simple  shear  of  DH-36  steel  with 
different  strain  rate;  2,5,11  -  tests  numbers,  b)  Stress  dependence  of  mean  temperature 
change  for  various  strain  rates  of  DH-36  steel;  temperature  does  not  vary  with  the  rate  of 

deformation  equal  10'  s'  . 

For  the  DH-36  steel,  temperature  does  not  change  in  elastic  range  of  shear,  i.e.  any 
thermoelastic  effect  was  observed.  The  rapid  increments  of  temperature  are  observed  in  the 
region  of  plastic  deformation.  The  temperature  increases  with  increase  of  the  rate  of 
deformation.  During  unloading,  the  temperature  smoothly  decreases  due  to  the  heat 
transmission  to  the  grips  and  surroundings. 

2.2.  Dynamic  tests  in  simple  shear  at  high  strain  rates.  The  device  with  double-shear 
specimen  is  placed  between  two  bars  of  the  SHPB  (Fig.  10).  The  mechanical  impedance  of 
the  shear  device  and  the  SHPB  must  be  the  same  to  avoid  the  noise  in  the  interface  signal. 
The  impulse  is  created  by  the  third  projectile  bar:  the  usual  compression  technique.  The 
impact  velocity  is  measured  by  the  set-up  consisting  of  three  sources  of  light  and  three 
independent  photodiodes.  The  time  interval  during  crossing  of  a  projectile  by  two  rays  of  light 
is  recorded  by  the  time  counter.  We  have  to  register  the  input,  the  transmitted  and  the 
reflected  impulses:  s, ,  st  and  sr  -  cf.  Fig.  11.  Impulses  are  registered  at  the  strain  resistance 
gauges. 


13 


Figure  10.  Standard  SHPB  apparatus.  P  -  projectile,  Bi  -  incident  bar,  Bt  -  transmitting  bar,  S 
-  shear  device  with  double  shear  specimen,  D  -  damper,  L  -  source  of  light, 

F  -  photodiodes,  TC  -  time  counter,  Tj ,  T2  and  Ts  -  strain  gauges,  A  -  amplifier, 

DO  -  digital  oscilloscope,  PC  -  microcomputer,  R-xy  recorder  or  graphic  printer,  THC  - 
thermovision  camera,  THV  -  oscilloscope  of  thermovision  set 


t[ns] 


Figure  11.  Incident,  reflected  and  transmitted  waves 

The  specimen  is  deformed  not  only  between  the  grips,  in  the  gauge  section.  The  part  under  the 
grips  is  partially  deformed  too.  The  transversal  strike  lines,  marked  on  the  specimen  before 
the  test,  become  after  deformation  the  curves  with  the  strike  sections  -  Fig.  12a. 
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Figure  12.  Displacement  between  external  and  internal  grips 


Relative  displacement  of  external  and  internal  parts  of  the  shear  device  Al  is  a  sum  of  two 
terms:  Al  =  Alg  +  2  Als,  where  Alg  corresponds  to  the  deformation  in  the  simple  shear  defined 
as  y=  Alg/a0  ,  and  Als  correspond  to  the  value  of  sliding  under  the  grips. 

At  the  first  step,  the  shear  device  together  with  the  universal  testing  machine  has  been  used 
under  quasi-static  loading  for  the  verification.  Measurements  of  the  displacement  or 
deformations  must  be  precise. 

It  was  assumed  that  in  the  simple  shear,  there  is  no  change  in  the  cross  section:  Se  =  l, 
d  =  const.  Then  the  mean  shear  stress  is  defined  by  <Jn  =  F/2  Se.  The  estimates  of  the  sliding 
value  under  the  grips  Als  can  be  obtained  in  the  quasi-static  test  of  loading  unloading.  After 
such  a  test,  the  permanent  deformation  of  the  specimen  can  be  measured  using  an  optical 
microscope.  Comparison  of  the  optical  measurement  with  those  obtained  by  displacement  at 
the  end  of  the  unloading  process  yields  the  value  of  Als.  This  procedure  must  be  repeated  for 
the  given  range  of  shear  deformations  y . 

Sliding  values  under  the  grips  (Als),  after  quasistatic  simple  shear  up  to  100% 
deformation  and  unloading,  for  test  with  different  strain  rate  are  shown  in  Table  1.  Three 
specimens  were  used  to  obtain  the  average  value  for  every  rate  of  deformation.  For  every 
specimen  three  lines  were  taken  to  obtain  this  value  (see  Fig.  13.) 


Table  1 


Strain  rate 

Als  (mm) 

10'3 

0.30 

10'1 

0.32 

W 

0.35 

15 


Shear  zones 


“Is! _ 1_ 


Jo 


Figure  13.  The  transversal  strike  lines  on  the  specimen  applied  to  measure  changes  of 
geometry  of  the  shear  zones  according  to  the  method  shown  in  the  Fig.  12. 

For  the  DH-36  steel  significant  part  of  deformation  take  place  out  of  the  shear  zones  -  under 
the  grips.  In  the  quasistatic  range  of  deformation,  the  contribution  of  sliding  under  grips 
increases  a  little  when  the  shear  deformation  increases. 

Average  value  of  the  deformation  range  out  of  shear  zone  was  Aao  =  0.52  mm  (see  Fig.  12). 
This  value  was  not  changed  with  rate  of  deformation. 

Measurement  of  dynamic  deformation  at  high  strain  rates  in  the  case  of  simple  shear 
of  the  metal  sheets  is  more  complicated  as  in  the  quasi-static  loading.  In  the  case  of  dynamic 
deformation  with  the  SHPB,  determination  of  the  relative  displacement  Al  of  the  external  and 
internal  parts  of  the  shear  device  is  not  a  quite  simple  matter.  It  is  assumed  that  both  bars  are 
elastic  and  displacements  8 f  t)  and  82(f)  of  two  faces  of  the  shear  device  can  be  determined 
from  the  measured  elastic  waves  et ,  st  and  sr\ 

(2.1)  8,  it)  =  cQ  j  [e,  (r)  (r)]  dr,  82  (t)  =  cjet  (r)dr, 

0  0 

where  c0  is  the  longitudinal  wave  speed  in  the  bars. 

The  mean  shear  strain  f*  corresponding  to  a  relative  displacement  of  two  parts  of 

(2.2)  r  =AL  =  3_[<?i  (o-^w] 

a0  ao 

Now  one  must  take  into  account  the  displacement  correction  due  to  the  sliding  of  the 
specimen  under  the  grips.  The  mean  shear  strain  y*(t)  can  be  expressed  as 

(2.3)  f(0  =  f*“— K 

a„ 

where  Als  is  defined,  as  a  function  of  shear  strain  y* .  The  mean  strain  rate  is  given  by  the  first 
derivative  f  =  dy/dt. 
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The  mean  shear  strain  <x12  can  be  calculated  from  the  relation 

(2.4)  5r12=[FI(0  +  F2(0]/4^f 

where  F  ft)  and  F2(t)  are  forces  acting  on  two  ends  of  the  shear  device,  they  can  be  expressed 
as  functions  of  eft),  er(t)  and  eft)  using  the  Hook's  law.  Also,  the  mean  shear  strain  is  given  in 
the  following  form 

(2.5)  <t,2  =-E  f L[£  (,)  +  (,)  +  e  (,)] 

4  A, 

where  E0  is  the  Young's  modulus  and  An  is  the  cross  section  area  of  Hopkinson  bars.  It  must 
be  remembered  that  all  those  relations  have  been  derived  using  the  elementary  wave  theory  of 
longitudinal  elastic  waves. 

We  obtain  directly  the  hardening  function  <Jn(y)  of  the  examined  specimen  in  simple  plane 
shear.  Using  this  method,  we  can  determine  the  surface  /  (crn,y,y)  =  0  if  we  assume  that  the 
time  of  transition  of  waves  trough  the  specimen  is  very  small  compared  to  the  total  time  of 
loading  by  the  first  pulse.  The  result  of  multiple  reflections  of  the  loading  waves  in  the 
specimen  is  the  uniform  distribution  of  strains  and  stresses  in  the  whole  gauge  section. 

The  an  (y)  curves  for  high  strain  rates  are  shown  in  Fig.  14,  for  the  DH-36  steel 
sheets.  In  this  experiment  we  use  the  projectile  of  the  length  150  mm,  strikes  with  the  impact 
velocity  from  10  to  12  m/s.  We  observe  very  high  level  of  stress  for  small  shear  deformations. 


v 


Figure  14.  Shear  stress  t=  <Jn~  shear  strain  y  curves  obtained  at  high  strain  rates  of  the 

3  1 

order  of  1.7  x  10  s'  for  DH-36  steel  sheets 
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2.3.  Metallurgical  observations.  In  the  Fig.  15  microstructure  in  DH-36  steel  before  testing 
is  shown.  Average  diameter  of  grain  was  a  10  pm.  Using  metallographic  microscope 
observation,  changes  in  the  microstructure  in  the  shear  zones  deformed  to  40%  after  dynamic 
simple  shear  are  not  visible. 


In  the  Fig.  16  we  can  observed  that  the  specimen  is  deformed  not  only  between  the  grips,  in 
the  gauge  section.  The  part  under  the  grips  is  partially  deformed  too.  The  transversal  strike 
lines,  provided  from  mechanical  polishing  of  the  specimen  before  the  test,  become  after 
deformation  the  curves  with  the  strike  sections. 


Figure  16.  Displacement  between  external  and  internal  grips  under  metallographic 
microscope.  Length  of  the  shear  zone  is  equal  to  3  mm. 
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2.4.  Dynamic  tests  in  simple  shear  at  very  high  strain  rates.  The  highest  strain  rate  in 
the  specimen  can  be  obtained  using  only  one  bar  of  the  SHPB  system.  We  use  the 
transmitting  bar  only  (tube)  and  the  shearing  device  is  placed  in  front  of  this  bar.  This 
technique  is  based  on  the  principle  given  by  Dharan  and  Hauser  [13].  In  this  new 
experimental  technique,  the  flat  projectile  strikes  directly  the  double  shear  specimen  placed  in 
the  device  on  the  front  of  the  transmitting  bar,  as  shown  in  Fig.  17.  In  this  technique,  the 
transmitting  bar  has  the  form  of  a  tube.  We  have  to  register  only  the  transmitted  impulse  st 
and  the  velocity  of  the  projectile.  The  set-up  with  three  light  beams  makes  it  possible  to 
determine  the  acceleration/deceleration  of  the  projectile  just  before  impact  and,  thus,  the  exact 
value  of  v0  can  be  determined. 

The  axial  displacement  Al  imposed  on  the  double  shear  specimen  by  a  projectile  impact  can 
be  found  as  follows:  Al(t)  =  8i(t)  -  S2(t),  where  Si(t)  is  the  axial  displacement  of  the  internal 
part  of  the  shearing  device  with  respect  to  the  external  part,  82(f)  is  the  displacement  of  the 
Hopkinson  tube,  given  by 


(2.6)  8l  (t)  =  vat  and  S2  ( t )  =  c0\st  (r)dr 

0 

where  va  is  the  impact  velocity  of  the  projectile,  c0  is  the  elastic  wave  speed  in  the  Hopkinson 
tube,  and  sft)  is  the  transmitted  signal  in  the  Hopkinson  tube,  measured  by  the  strain  gauge 
Tj.  The  axial  displacement  8i(t)  of  the  central  part  of  the  double  shear  specimen  is  measured 
as  a  function  of  time  by  an  optical  extensometer,  acting  as  a  contact-less  displacement  gauge. 
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Figure  1 7.  Configuration  of  experimental  set-up  for  direct  impact.  S  -  shear  device  with  a 
double  shear  specimen,  HT  -  Hopkinson  tube  (transmitter  bar),  P  -  projectile, 

D  -  damper,  L  -  laser,  F  -  photodiodes,  TC  -  time  counter,  77  -  strain  gauge, 

A  -  amplifier,  DO  -  digital  oscilloscope,  PC  -  computer 


The  mean  shear  strain  and  the  mean  strain  rate  can  be  found  as 

1 


(2.7) 


m=— 


Vot~Co\£,(T)dT 


r(t)  =—[v0-c0£l(t)\ 
a„ 


The  shear  stress  can  be  expressed  as 


(2.8) 


ZoE 
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where  E0  is  the  Young  modulus  of  the  tube,  ST  =  n{d\  -df)/  4,  di  and  d?  are  the  internal  and 
external  diameters  of  the  tube,  respectively.  Thus,  the  shear  stress  in  the  double-shear 
specimen  is  proportional  to  the  current  signal  of  the  transmitted  longitudinal  wave  st. 

A  typical  oscillogram  obtained  from  the  direct  impact  technique  for  DH-36  sheet  metals  is 
shown  in  Fig.  18.  Because  of  a  high  strain  rate  of  the  order  y  =  6.2  -103  s~l ,  vibrations  in  the 
shear  device  are  inevitable. 


100 

0 

-100 

-200 

-300 

-400 

-500 

-600 

-700 

-800 


. _ i _ , 

J  SC202  45 1 

‘A . 

t 

... 

f 

. i 

/ 

y 

.  i  1 . 

/ 

\  a/S 

i 

i 

200 
t  [\xs] 


Figure  18.  Oscillogram  from  the  direct  impact  test  at  velocity  vQ=  20  m/s 


Results  of  experimental  investigations  of  DH-36  steel  in  the  case  of  dynamic  simple 

3-1  3-1 

shear  at  very  high  strain  rates  from  5.2  10  s'  to  6.2  10  s'  are  presented  in  the  Fig.  19.  We 
can  compare  shear  stress  -  shear  strain  curves  obtained  in  the  dynamic  test  at  very  high  strain 
rates.  From  this  picture  we  are  able  to  determine  the  shear  strain  corresponding  to  ductile 
fracture  from  the  condition  dx/dy  =  0. 


Figure  19.  Stress-strain  relation  obtained  during  dynamic  simple  shear  of  DH-36  steel 
(direct  impact);  189  202  -  tests  numbers 
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Figure  20  shown  the  comparison  of  the  shear  stress  -  shear  strain  curves  in  case  of  the 

'J  1 

dynamic  loading  (direct  impact  test)  at  the  strain  rate  6.3  10  s'  with  the  quasi-static  curves. 


Figure  20.  Comparison  of  stress-strain  characteristics  obtained  during  quasistatic  and 

dynamic  tests 


In  the  case  of  quasi-static  simple  shear  of  DH-36  steel  the  maximum  values  of  shear 
stresses  are  situates  on  the  level  of  400  MPa  -  cf.  Figs  7  to  9.  The  comparison  of  maximum 
shear  stress  in  the  quasi-static  and  dynamic  tests  at  high  strain  rate  higher  then  103  s'1  is  show 
on  the  Fig.  21. 


log(dy/dt) 

Figure  21.  Shear  stress  r  versus  logarithmic  shear  strain  rate  dy/dt 
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In  the  Table  2  all  experimental  results  in  dynamic  test  (in  the  case  of  SHB  -  Split  Hopkinson 
Bar  technique  and  in  the  case  of  DI  -  direct  impact  of  the  projectile  into  the  shear  sample)  are 
collected. 


Remarks  to  the  Table  2. 


scale 

Generally 

shear  bands 

0 

lack 

absent 

1 

Weak 

matt  at  the  specimen  surface 

2 

Strong 

sharp  line  of  the  shear  bands  on  the  specimen  surface 

Impactor 

Length 

LP 

L=300  mm 

MP 

L=150 

Test 

DI  direct  impact 

SHB  Shear  in  Split  Hopkinson  Pressure  Bars 
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Table  2.  Results  of  dynamic  tests  in  simple  shear 


Test 

number 

Specime 

n 

number 

Impactor 

Test 

a0  [mm] 

vp  [m/s] 

Ekp  [J] 

21 

SHB 

0.650 

5.45 

4.4 

16 

SHB 

0.640 

0.0 

SHB 

0.647 

6.94 

7.2 

SHB 

3.03 

1.4 

24 

SHB 

0.624 

10.65 

17.0 

164 

18 

MP 

SHB 

0.635 

10.50 

16.5 

19 

SHB 

0.640 

10.04 

15.1 

22 

SHB 

0.654 

25 

SHB 

0.656 

26 

SHB 

0.652 

12.88 

24.9 

27 

SHB 

0.640 

11.98 

21.5 

170 

28 

SHB 

0.648 

3.97 

2.4 

29 

SHB 

0.650 

13.96 

29.2 

30 

SHB 

0.651 

14.63 

32.1 

SHB 

0.642 

13.88 

28.9 

SHB 

0.651 

177 

LP 

SHB 

0.656 

7.96 

19.0 

34 

LP 

SHB 

0.655 

8.03 

19.4 

LP 

SHB 

0.650 

7.57 

17.2 

17 

Dl 

0.635 

16.99 

43.3 

Dl 

0.647 

16.37 

40.2 

37 

Dl 

0.648 

19.69 

58.2 

Dl 

0.650 

19.32 

56.0 

40 

Dl 

0.652 

19.67 

58.1 

38 

MP 

Dl 

0.630 

3.77 

2.1 

Classification  trees  are  used  to  predict  membership  of  cases  in  the  classes  of  categorical 
dependent  variable.  To  calculating  we  have  used  application  the  STATISTICA  version  6.  All 
results  of  dynamical  experiments  are  gathered  in  Table  2.  As  dependent  variable  the  partial 
cracking  was  chosen.  For  physical  reason  we  have  get  a  set  of  independent  variables  listened 
in  the  Table  3. 


Table  3.  Importance  of  the  independent  variables 


Variable 

rank 

Importance 

kinetic  energy  of 
impactor 

Ekp  [J] 

Ekp  [J] 

82 

0.817353 

deformation 

Tend 

gend 

63 

0.626555 

deformation  rate 

dy/dt  max 

dg/dt  max 

100 

1.000000 

maximal  shear 
stress 

Tmax  [MPa] 

tmax  [MPa] 

100 

1.000000 

average  shear 
stress 

^avg  [MPa] 

tavg  [MPa] 

64 

0.642952 

100  -  the  best 

Variable  rank  is  calculated  by  algorithm  of  the  classification  using  correlation  between  all 
variables.  The  classification  algorithm  was  applied  to  N=23  samples.  According  this 
classification  the  most  important,  strong  correlated,  are:  deformation  rate,  maximum  shear 
stress  and  weaker  kinetic  energy  of  projectile,  see  Table  3  and  Fig.  22.  In  the  Fig.  23  levels  of 
the  classification  for  both  the  most  important  variables  are  presented,  dy/dt=5800  1/s  and 
'tmax=369.5  MPa.  Algorithm  of  classification  is  stopped  on  the  second  level  of  the 
classification  tree  when  the  split  error  of  the  samples  is  less  then  0.05. 


Figure  22.  Importance  plot.  Dependence  variable:  partial  cracking 


- 0 

—  1 


Figure  23.  Tree:  graph  for  variable  partial  cracking.  Number  of  non-terminal  nodes  2, 

Number  of  terminal  nodes  3 


2.5.  Comparison  of  simple  shear  test  with  compression  test.  To  understand  the  mechanical 
response  of  DH-36  steel,  uniaxial  compression  tests  are  performed  on  small  cylindrical 
specimens,  using  an  Instron  testing  machine.  True  strains  exceeding  40%  are  achieved  in 
these  tests  in  the  quasistatic  conditions  of  loadings,  with  strain  rates  of  the  order  of  6  10'  s'  . 
All  tests  are  carried  out  on  a  DH-36  baseplate.  All  samples  used  in  the  compression  test  have 
8  mm  nominal  diameter  and  8  mm  heigh.  To  reduce  the  end  friction  on  the  samples  during 
deformation,  the  sample  ends  are  polished  and  greased.  The  stress-strain  relation  in  the 
compression  test  is  presented  in  the  Fig.  24. 


-2  -1 

Figure  24.  Quasi-static  compression  test  of  the  DH-36  steel,  strain  rate  5.77  10'  s' 
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Comparison  of  the  result  obtained  in  the  simple  shear  test  and  the  compression  test  for 

2  1 

DH-36  steel  at  low  strain  rate  5.77  10'  s'  is  shown  in  Fig.  25.  We  assume  the  Huber-Mises- 
Hencky  yield  criterion:  y  =  £u  /a/3  and  <x12  =-j3ern  . 

In  fact,  the  influence  of  stress  components  ai  /  and  U22  must  be  taken  into  account.  The  stress 
intensity  is  defined  as:  cr  =  [cr,2,  +  <r22  +  3a 21  -  crna22  ]1/2 . 

We  can  observe  that  the  hardening  of  the  material  is  similar  for  different  values  of  the  strain 
rate  (the  same  effect  is  observed  in  the  case  of  the  Tresca  yield  condition).  In  the  case  of 
Huber-Mises-Hencky  (HMH)  yield  condition,  an  important  difference  in  stress  increment  is 
observed.  The  existence  of  a  difference  between  simple  shear  tests  and  tensile  tests  after 
application  of  the  HMH  criterion  was  also  observed  by  various  authors.  They  explain  this 
difference  by  the  different  systems  of  slip  activated  according  to  the  type  of  loading  and 
interaction  of  displacements  at  long  and  short  distance.  Also,  for  the  metal  like  the  DH-36 
steel,  the  HMH  yield  criterion  is  not  the  best. 


Figure  25.  Comparison  of  the  shear  and  tension  test  at  strain  rate  5.77  10 
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3.  Constitutive  model  in  application  for  finite  element  simulation  of  dynamic 
experiments  and  construction  of  the  dynamic  fracture  locus 

The  simple  shear  in  the  direction  ei  of  the  coordinate  system  (e  i ,  ei)  is  defined  by  the  relations 
u\ =  y(f)  x2  ’  u2  =  u3  =  o , 

Vi  =yx2  ,  v2=v3=0, 

where  y  =  tg  cp  (cf.  Fig.  26)  and  y  are  the  plastic  shear  strain  and  shear  strain  rate, 
respectively. 

*2 


Figure  26.  Scheme  of  simple  shear 


From  the  velocity  field  v,  the  velocity  gradient  V,  the  rate  of  deformation  D  and  the  material 
spin  (o  can  be  determined  in  the  system  (ei,  ei)  as 


"0  f 

y  0  f| 
D  =  — 

y  0  1) 

CD  =  — 

vO  0, 

2U  OJ 

2-1  OJ 

The  Cauchy  stress  tensor  a  has  the  non-zero  components 


The  program  of  finite  element  method  was  used  to  the  numerical  simulations  of  the 
formulated  problem  of  quasi-static  and  dynamic  simple  shear  of  thin  sheets.  We  assume  the 
initial  and  boundary  conditions  similar  to  those  used  in  the  experiment. 

In  the  case  of  quasi-static  loading,  the  boundary  conditions  have  the  following  form: 
in  fixed  ends  of  specimen  (under  grips)  -  cf.  Figure  26: 

ufxl,x2,t)  X2=0  =  u2(x1,x2,t)  Xi=0  =  o, 
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where  vQ  is  the  velocity  of  the  testing  machine  in  traction.  In  this  case,  the  sliding  of  material 
under  the  grips  is  neglected. 

At  the  free  ends  of  the  specimen  (for  x\=  0  and  xi  =  la )  we  have: 


=  0. 


x,  =0 

*i=t> 


In  the  case  of  dynamic  deformations,  these  conditions  must  be  replaced  by  conditions  of 
balance  of  forces  in  the  contact  between  the  specimen  and  the  measuring  bars.  We  should 
remember  that  the  shear  device  and  bars  of  SHPB  system  have  identical  mechanical 
impedances. 

In  the  analysis,  we  must  take  into  account  that  the  loading  of  the  specimen  is  not 
instantaneous.  The  loading  compression  wave  must  take  some  time  to  transmit  from  one  end 
of  the  device  to  the  other.  In  the  first  steep  a  simplifying  assumption  can  be  introduced,  and 
the  process  of  waves  propagation  in  the  specimen  can  be  neglected.  However,  it  is  worthwhile 
to  notice  that,  in  our  tests,  a  very  accurate  equilibrium  is  obtained  regarding  the  forces  on  two 
sides  of  the  shear  device,  in  the  time  period  50  ps  <  t  <  500  ps.  The  force  is  taken  to  be  equal 
to  the  mean  value  of  input  and  output  force. 

The  amplitude  of  loading  is  determined  from  the  dynamic  experiment.  We  assume  that  in  the 
contact  between  the  specimen  and  measuring  bars,  the  force  is  constant  in  time  and  equal  to 
Fmax.  The  boundary  conditions  take  the  form: 


cr12 (*j ,x2,t)  X2=0  =  -Fmax  /  SE  and  u2 (x, ,  x2 ,  f)  ,2=0  =  0 . 


max 


We  assume  the  homogeneous  zero  initial  conditions. 

In  the  finite  element  method  the  rectangular  mesh  is  introduced.  Deformation  of  the  mesh  in 
time  is  determined.  At  the  same  time,  the  components  of  the  stress  tensor  an  and  a>n  = 


the  stress  intensity  at  =  (3/2  Sy  Sy)m  and  the  equivalent  strain  e,  =  (2/3  s/  Sify  are 
determined. 
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3.1.  Constitutive  description  of  inelastic  deformation  in  application  for  the  study  of 
strain  rate  effect  on  ductile  fracture.  We  assume  that  a  multiscale  approach,  which  was 
proposed  in  relation  to  plastic  flow  of  metals  accounting  for  shear  banding,  in  refs  [17-19], 
can  be  extended  to  describe  plastic  deformation  of  DH-36  steel.  This  is  justified  by  the  results 
of  microstructural  observations  presented  in  [15],  where:  “shear  bands  are  seen  to  have 
formed  at  an  angle  45°  to  the  loading  direction  (Figs.  17  and  18),  and  within  the  shear  bands, 
the  grains  are  seen  to  have  been  extended  in  the  shear  band  direction  ”,  p.  1034  in  [15], 

Among  many  possible  realizations  of  shear  banding  processes  one  can  single  out  the 
group  of  processes  characterizing,  at  lest  approximately  or  for  sufficiently  long  deformation 
paths,  with  the  same  contribution  of  two  systems  of  shear  bands,  i.e.  ff  =  fff ,  where  f'SB  is 

the  shear  banding  function  describing  the  contribution  of  the  shear  banding  system  /Wand 
i=2,  respectively,  in  the  total  plastic  shear  strain  rate  [19].  On  the  other  hand,  the  non- 
symmetric  activation  of  shear  bands  can  be  induced  by  a  sufficiently  large  change  of  loading 
direction  or  rotation  of  the  principal  axes  of  the  stress  tensor  in  the  inhomogeneous 
deformation  process  caused  by  boundary  conditions.  Therefore,  we  postulated  in  [20]  that: 


(3-1)  A fSB  =0  for  8  e  [0,  8C]  and  A fSB  =  zl(S)  for  8 


x  71 
2 


,T':pF 


cos  8. 


The  function  A(b)  being  a  measure  of  the  mentioned  asymmetry  of  shear  banding  should  be 
identified  from  numerical  simulations  of  the  experiment  accounting  for  the  change  of  loading 
direction. 

Because  of  the  above  assumptions,  the  plastic  flow  laws  take  the  form: 


•  For  the  case  in  which  the  loading  direction  described  by  the  objective  rate  of  stress  t 

r 


is  pointing  at  partially  active  range,  i.e.  for  8 ' 


x  n 
S‘’2 


(3.2) 


1)  = 


t  :  iiF 

- Llr  + 

m-fsB) 


x:pF 
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)  tan  20 
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•  For  the  case  in  which  the  objective  rate  of  stress  t  is  pointing  at  fully  active  range,  i.e. 
for  8  e  [0,8C]: 


(3.3) 


I)  = _ _ 

m-fsi) 


M/  = 


42k 


7  1 

The  Huber-Mises  yield  condition  is  assumed:  J2  =  k  ,  where  J2  =  —  t':  t'  and  the  symbol  k 

denotes  the  yield  limit  obtained  in  a  shear  test.  The  above-mentioned  assumptions  leading  to 
the  simplified  formula  (3.3)  find  confirmation  in  experimental  observations  for  metallic 
materials,  which  show  that  the  spatial  pattern  of  micro-shear  bands  does  not  change  for 
loading  conditions  that  deviate  within  limits  from  the  proportional  loading  path,  i.e.  the  load 
increments  are  confined  to  a  certain  cone  (fully  active  range)  the  angle  of  which  can  be 
determined  experimentally,  [19].  The  specification  of  these  equations  and  the  identification  of 
the  shear  banding  contribution  function  fSB  are  discussed  in  [20]. 

If  the  dependence  on  strain  rate  comes  into  play,  the  viscoplasticity  flow  law  can  be 
applied  [21]: 


(3.4)  T>p  =  yvppF ,  pF  =  - -  t'  , 

where  yvp  is  the  viscoplastic  shear  strain  rate  determined  in  the  following  way  for  the 
exponential  overstress  function 


(3.5) 


■vp 


=  Y  o  exp 


h 


Where  the  power-like  strain  hardening  law  is  assumed. 

According  to  the  analysis  in  [20],  the  shear  banding  function  takes  the  form: 


(3.6) 


f  SB  ~ 


fsBn 


1  +  exp(«  —b  yvp) 
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The  yield  condition  J2  =  kd  is  obtained  by  inverting  (3.5) 


(3.7)  kd  ={l-fSB\A  +  B(fp  )n\l  +  Clni*\ 


yvp 

where  y*  =- — ,  y0  =2xlOw s~! . 

Y  o 

The  formula  (3.7)  shows  that  Perzyna  model  with  the  exponential  overstress  function  and  the 
assumed  strain  hardening  law  is  equivalent  to  the  Johnson-Cook  model  considered  in  [15]. 

In  our  case  the  results  presented  in  [15]  can  be  used,  where  the  Johnson-Cook  model 
was  identified  for  DH-36  steel.  If  the  effects  of  shear  bands  and  temperature  are  neglected  we 
have 


(3.8) 


(A  +  B(y^)")(l  +  Clni’), 


where  A  =  1020MPa, 


B  =  1530, 


C  =  0.015, 


*  yvp 

0.4  and  y  =  - — ,  y0 
Y  o 


=  2x10 


10  -i 


3.2.  Conclusions 

Application  of  the  general  scheme  of  calibration  of  fracture  model  proposed  e.g.  by  Y.  Bao, 
Y.-W.  Lee  and  T.  Wierzbicki  [16]  to  the  investigated  experimentally  material  (DH-36  steel) 
consists  of  the  following  steps: 

1.  Experimental  determination  of  stress-strain  characteristics  of  DH3-steel  under  quasi-static 
and  dynamic  loading  conditions,  e.g.  in  tension  or  compression  tests. 

2.  Identification  of  viscoplasticity  model  accounting  for  shear  banding. 

3.  Measurements  of  force-displacement  response  of  simple  shear  test. 

4.  Development  of  finite  element  model  of  a  specimen  and  numerical  simulation  of  the 
loading  process. 
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5.  Calculation  of  the  evolution  of  equivalent  strain  and  the  stress  triaxiality  at  the  critical 
point  (detected  during  the  test  with  use  of  thermovision  camera). 

6.  Determination  of  the  equivalent  strain  to  fracture  from  the  comparison  between  the  tests 
and  numerical  simulations. 

An  example  of  the  force-displacement  response  of  simple  shear  test  performed  with  use  of 
DH-366  steel  samples  is  shown  in  Fig.  27.  The  development  of  finite  element  model  of  a 
specimen  used  in  the  simple  shear  test  is  depicted  in  Fig.  28.  This  makes  a  basis  for  the 
numerical  simulation  of  the  loading  processes  at  different  strain  rates  and  performing  of 
further  steps  of  the  calibration  procedure.  Realization  of  these  steps  should  contribute  in 
future  to  the  construction  of  the  limiting  failure  curve  (fracture  locus)  for  different  strain  rates. 


Figure  27.  Load  -  extension  curves  obtained  in  simple  shear  tests  for  different  strain  rates, 

DH-36  steel. 
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Figure  28.  Assumed  shape  of  the  specimen  and  the  initial  mesh. 


4.  Annex  I  -  PP  -  presentation  -  seminar  at  MIT 

Prof.  R.  B.  Pqcherski  took  part  in  the  Workshop  on  Computational  Fracture  Mechanics  and 
ICL  Review  Meeting,  MIT,  October  13-14,  2005  and  presented  the  lecture  prepared  jointly 
with  Prof.  W.K.  Nowacki  on:  “A  New  Experimental  Technique  for  Studying  Strain  Rate 
Effects  on  Fracture”.  During  the  lecture  the  potential  in  the  field  of  experimental  and 
theoretical  investigations  of  the  team  at  the  Institute  of  Fundamental  Technological  Research 
in  the  joint  research  project  was  presented.  The  main  subject  of  the  discussion  after  the 
presentation  was  related  with  the  simple  shear  test  at  different  strain  rates  and  the  possibilities 
of  detecting  the  region  of  fracture  initiation. 

5.  Annex  II  -  Micromechanics  of  localized  fracture  phenomena  in  inelastic  solids  -  by  Z. 

Nowak,  P.  Perzyna  and  R.B.  Pqcherski 
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Our  tradition  in  fracture  research 

Phenomenological  criteria  of  strength  (material  effort): 

-  energy  based  approach  for  isotropic  solids  (M.T.  HUBER 
[1904] 

-  Volume-distortional  energy  hypothesis  accounting  for  the 
effects  of  triaxial  states  -  W.T.  BURZYNSKI  [1928], 

-  application  of  energy  based  approach  for  orthotropic 
solids  -  W.T.  BURZYNSKI  [1928], 

-  Energy-based  approach  for  solids  with  arbitrary  symmetry 
J.  RYCHLEWSKI  [1984], 

-  Energy-based  apprach  applied  for  cellular  solids  and 
foams  of  different  symmetry  R.B.  PECHERSKI  et  al. 
[2003-2005]. 


Industrial  Fracture  Consortium  Meeting,  MIT  Cambridge,  October  13-14,  2005 
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M.  T.  HUBER  (1872-1950) 
Czasopismo  Techniczne, 
Lwow  (Lemberg)  [1904] 


Bys.  I.  Kozklad  napr^zen  rv  ciele  odksztaleouym:  a)  stall  naturalny,  b)  roze 

w  kierunku  y,  c)  rozci^ganie  w  kierunkacb  x  i  y 


Prwfw»ir  M  T  I  lie  Rrrtm  if  L.wif  PalYtMfcmr  .IMll 


’’Specific  Strain  Work  as  a  Measure  of  Material  Effort  -  A  Contribution 
to  Foundations  of  the  Theory  of  Material  Strength’’-  ARCHIVES  OF 
MECHANICS,  56,  No.  3,  173-190,  2004  (English  translation  on  the 
centenary  of  the  first  Polish  edition). 
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The  specific  work  of  strain  is  equal  to  the  work 
performed  by  intermolecular  forces  acting  between 
molecules  ....  Since  the  work  of  intermolecular  forces  is 
the  greater  the  larger  are  the  relative  displacements  of  the 
molecules  in  strain  a  very  probable  hypothesis  suggests 
itself  that: 

Material  effort  is  measured  by  specific 

work  of  strain. 

Then,  if  in  a  certain  point  of  the  body  work  of  strain  goes 
beyond  a  determined  value  depending  on  the  material,  a 

permanent  separation  of  molecules  of  the  body,  that  is 
its  fracture,  must  occur. 


Industrial  Fracture  Consortium  Meeting,  MIT  Cambridge,  October  13-14,  2005 
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Huber-Mises  criterion 

Linear  elasticity  -  Hooke’s  material: 
strain  work  =  stored  elastic  energy 


<D  =  Ov  +  <3>f 


R.von  Mises  (1913)  -  quadratic 
approximation  of  Tresca  criterion. 

Richard  von  MISES  (1883-1953) 


Industrial  Fracture  Consortium  Meeting,  MIT  Cambridge,  October  13-14,  2005 


6 


W.  T.  BURZYNSKI  (1900-1970) 

„Study  on  Material  Effort  Hypotheses”, 

Lwow  1928. 

„Ueber  die  Anstrengungshypothesen”, 
Schweizerische  Bauzeitung,  94,  259-162,1929 

d>  y  -h  r)  O  v  =  O  gr  =  const ,  0  <r\<  1 
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2 
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co2+r| - co7  =  const 


-  v2  y + (^2  -  ^5  y + (<*i  -  ) 
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W.T.  Burzynski  [1928]  -  R.  Hill  [1949] 
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Energy  based  approach  for  orthotropic 
solids  -  W.  Burzynski  [1928] 

0  =  0,,  +  O/ 

cDv(a)  =  \B*(o  11  +  022  +  033)2 

0/(a)  =  J-[7V*(0U  -  a22)2  +Z,*(a22  -  a33)2  +  M*(033  -  an)2] 

+  2  (i5*023  +  g*C^3  +  i?* 0j2  ) , 

5*  =  c„  +  c22  +  c33, 


L*  =  \{B*  -C33), 
M *  =  2-C8*  -C32), 
AT  =  \{B*  -Cl2), 


4 P*  =  C  44  +  2C45 


4g*  =  C55  +  2C56#^ 

<--26 


4i?*  =  C66  +  2C64 


C 16 
C34  ' 


M.  Janus-Michalska  &  R.B.  Pecherski 
Technische  Mechanik,  23,  2003. 

Criterion  for  metal  foams:  C.San  Marchi 
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Division  of  Applied  Plasticity 

Head:  Prof.  Wojciech  K.  NOWACKI 

v 

Thermo-mechanical  modelling: 


•  constitutive  description, 

•  numerical  identification  &  verification 
Materials:  steels,  Dual  Phase  (DP)  steel,  TRIP  steel, 
intermetallics 
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Experimental  investigations: 

•  Simple  shear  tests  at  different  strain  rates  from 

•  10"3  s'1  to  104  s'1 

•  Compression  and  tension  tests  at  different 
strain  rates  from  10'3  s1  to  103  s1 

•  T  aylor  test 

•  Kid(t),  crack  propagation  velocity 

•  Damage  in  brittle  materials 

•  Ductile  fracture 

Materials:  steels,  polymers,  SMA,  composites 
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Testing  machine  25  kN  and  the  compression  split  Hopkinson  bar  facility 

Industrial  Fracture  Consortium  Meeting,  MIT  Cambridge,  October  13-14,  2005  13 
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I.  Compression  Split-  Hopkinson  Pressure  Bar  (SHPB)  scheme 

dynamic  experiment:  a)  simple  shear 

b)  tension 


Strain  gauge 


Strain  gauge 


input  bar  shear  device  output  bar  damper 
(tension  device) 

IMPORTANT:  The  mechanical  impedance  of  two  parts  of  the  shear  device 
and  SHPB  are  the  same. 


s  ,  £  ,£  and  v 

i  '  r  '  t  o 
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II.  Schematic  of  configuration  of  experimental  set-up  for  direct  impact: 

a)  simple  shear  test  b)  tension  test 


We  use  the  transmitted  bar  only  and  the  shearing  (tension)  device  is  placed 
in  the  front  of  this  bar. 

-pj  'aser  strain  gauge 


photodiodes 


shear  device 
(traction  device) 


output  bar 


damper 


W 


Pd  ’  Ad  <=>  Pshpb  ’  4m 
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In  the  static  investigations  of  the  thin  sheets  we  can  distinguish  two 
approaches,  based  on  different  principles  of  the  shear  device. 


a)  Pure  shear 


rigid  specimen  clamp 
d. 


rigid  specimen 
'  clamp  with  a 
horizontal  and 
vertical  carriage 


specimen 


t  shear  zone 


b)  Simple  shear 


symmetry  axis. 


external 

part 

d 


\ 


-internal  part 


I  ^2F 


rix' 


external 
’  part 


specimen 


t  shear  zones t 


Shear  device  of  the  one  shear  zone  only 
proposed  by  G’Sell  (1983) 


Shear  device  for  ..double  shear”  specimens 
proposed  by  Myauchi  (1978) 
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shearzones  a  =  3  mm 


double  shear  specimen 
(thickness  of  d=  0.8  mm) 


The  greatest  double  shear  specimen  height 
allowed  by  the  shear  device  is  lQ  =  30  mm. 


The  width  value  a0  is  chosen  with  a  view  to 
satisfying  two  requirements: 

•  the  ratio  aJ6  of  the  width  to  the  thickness  has 
to  be  small  enough  to  avoid  the  buckling  of  the 
sheared  zone  ( aJ6<  10) 

•  the  ratio  aj\0  of  the  width  to  the  height  has  to 
be  sufficiently  small  to  minimize  the  error  due  to 
the  non  homogeneity  of  the  shear  stress  and 
strain  at  the  two  ends  of  the  sample 

(a  J!  <1/10). 
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Simple  shear  -  controlled  by  a  SHPB  acting  in  compression 


time  t  =  0  h0 


_ ■ 

incident  bar 


time  t1  >0 


transmitter  bar 


incident  bar 


time  t2  >  time  t1 


Device  for  high  strain  rate 
shear  test 
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A  compression  split  Hopkinson  pressure  bar  facility  at  Thermoplasticity  Laboratory 
(IPPT-  Warsaw) 


- ►  Bars  of  the  diameter: 

•  22  mm  -  steel 

•  25mm  steel  and  40mm  PA7 

- ►  Tube  of  the  diameter  40  mm 

Diameter  of  the  projectile  22  mm 


•Simple  shear  tests 


•Compression  and  tension 
tests 


•Taylor  test 


•Kid(t),  crack  propagation 
velocity 


V^Damage  in  brittle  materials. 
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A  new  compression  split  Hopkinson  pressure  bar 


Bars  of  the  diameter: 

•  25  mm  -  steel 

•  25mm  steel  and  40mm  PA7 
Tube  of  the  diameter  40  mm 
and  60  mm 

Diameter  of  the  projectile  25  mm 
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Simple  shear  in  the  classic  configuration  of  compression  SHPB 


\ 

The  input  force  and  the  output 

are  very  similar  in  shape. 

_ _ _ ! _ / 
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Configuration  of  experimental  set-up  for  direct  impact  -  Hopkinson  tube 


Experimental  stand  at  IPPT  (Warsaw) 


L 


I 

I 

I 

I 

I 

I 

I 

I 


I 

I 

I 


I 


The  temperature  of  the  surface  of  the  specimen 
can  be  measured  only  along  one  line. 


Vot~Co\£t  (T)dT 
0 


r(t)=  k  c0s,(t)\ 

vl2U)  =  ^rEo6)(t) 

ao 

2  SE 
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One  of  shear  paths  after  dynamic  simple  shear  at  strain  rate  3x1 03  s-1. 


A  photograph  of  the  circular  marking  before 
(diameter  of  circles  1  mm) 

and  elliptic  after  plastic  deformation  in 
simple  shear  up  to  0.91 . 


Within  the  central  part  of  the  shear  zone  3.0mm 
the  deformation  field  was  quite  uniform. 

Free  and  of  the  specimen  in  simple  shear. 

One  can  observe  the  heterogeneity  of  the 
displacement  field  and  important  sliding  under 
grips. 
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The  layers  under  the  grips  are  also  partially 
deformed.  The  transversal  straight  lines 
marked  on  the  specimen  before  the  test, 
after  deformation  show  small  curved 
transition  zones  with  the  straight  segment 
in  the  middle. 


Relative  axial  displacement  of  the  external  and  internal  grips  of  the  shear  device  is  relatively  easy 
to  measure  during  the  quasi-static  deformation  in  the  testing  machine.  It  is  a  sum  of  two  terms: 

Al  =  Alg+  Als 

where  Als  value  of  sliding  under  the  grips. 
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Macro-  and  micro-bands  of  shear  in  the  specimen 

An  essential  feature  is  the  formation  of  the  tangled  structure  and  of  the  dislocations 
cells.  Their  elongation  and  arrangement  tend  to  be  aligned  along  the  shear  direction. 


At  the  very  high  strain  rate  the  shear  macro-bands  are  observed. 


We  can  suppose  that  it  is  a  critical  strain  at  which  the  shear  localization  occurs.  Before 
arriving  at  the  critical  strain,  the  deformation  is  homogeneous  over  the  whole  gauge 
length  of  the  specimen. 


direction 

of 

traction 
C  x  axis  ) 


Shear  deformation  of  the  order  of  73%. 

In  several  grains  the  micro-bands  of 
shear  parallel  to  the  direction  of  x1  axis 
are  observed. 


ES  steel:  C  -50,  Ni  -25,  Cr-18,  Mn  -189,  Cu  -23,  Al  -57,  Si  -4,  P-17  in  10‘3  volume  percent 
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Shear  stress  -  shear  strain  relations 
strain  rates  from  10-3  s_1  to  102  s_1 
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Comparison  of  the  simple  shear  test  with  the  simple  compression  test  for  ES  steel  sheets. 


Maximal  stress  (MPa) 
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□  quasi-static  and  dynamic  test  in  compression 

■  simple  shear:  Gary  and  Nowacki  (1994) 

+  shear  test  in  hydraulic  machine 

•  direct  impact:  Klepaczko,  Nguyen  and  Nowacki  (1999) 


Variation  of  the  maximal  stress  versus 
logarithmic  rate  of  deformation,  for  different 
kinds  of  experiments. 


Rate  sensitivity  of  the  maximal  stress  for 
ES  steel  versus  logarithm  of  shear-strain 
rate. 
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Hopkinson  bar 
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FORCE 


TEMPERATURE 


r 


THERMOVISION 

CAMERA 


ELONGATION 


TIME 


TESTING  MACHINE 


Thermographic  set¬ 
up  to  measure  the 
temperature 
increase  with 
plastic  deformation 
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Theoretical  simple  shear  analysis 

The  simple  shear  in  the  direction  e.,  of  the  co-ordinate  system 
(e1  e2)  is  defined  by  the  relations 

u,=r(t)x2,  m,  =  m3  =  o,  \=yx2,  v2=v3=o, 

where  y-  tg  cp-  plastic  shear  strain,  /  -  shear  strain  rate. 


y  ( 0 
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T  =  /?LD 


3^D-(T-n) 

(72y  H 


[(t  -  n) + p] 


j 


1  if  /  =  0  and  D-(T-II)>0, 

0  if  /  =  0  and  D(T-II)<0  or  /<  0 


where  T  =  J3g,  j3  =  p/p ,  T  =  T  -  coT  +  Tco 
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Huber-Mises-Hencky  yield  criterion: 

Where  cry  is  the  yield  stress  in  simple 
a  corresponds  to  the  size  of  the  yield  surface:  a  =(T-n)  Dp . 


tension,  &  is  the  temperature  and 


The  change  of  temperature  is  described  as: 


The  term  on  the  right  hand  side  represents  the  rate  of  energy  dissipation  and, 
therefore,  n<\.  For  numerous  metals  stakes  the  value  from  0.02  to  0.1. 

30.00  — | 

For  plasticity  with 

20.00  —  /  k 

kinematic  hardening 

H  =l+— 

2// 
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Numerical  simulation  of  the  experiment  -  quasi-static  deformation  at  10-3  s 1 
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Adiabatic  shear  bands  -  ES  steel  (French  Standards) 
strain  rate  -  10-3  s_1 


Results  of  numerical 
simulation. 

(1/2  of  the  specimen 
in  simple  shear) 


Experiment 
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Numerical  simulation  of  the  experiment  -  dynamic  deformation  at  103  s 1 


Equivalent  strain 


y=5.1%, 


7=9.0%, 


7=18%, 


7=27% 


7=52% 


7=79% 
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Numerical  simulation  of  the  experiment  -  dynamic  deformation  at  103  s 1 


y=5.1%, 


y=9.0%, 


y=18%, 


Equivalent  stress 


Industrial  Fracture  Consortium  Meeting,  MIT  Cambridge,  October  13-14,  2005 


34 


Institute  of  Fundamental  Technological  Research 

f - N 

Conclusions  from  numerical  simulation  of  the  experiment 

v _ J 


Quasi-static  case: 


One  can  observe  the  heterogeneity  of  the  strain  and  stress  fields  at  the  ends  of 
the  specimen,  at  the  distance  less  then  1 .7%  (accurate  to  ICH  of  deformation)  of 
the  total  length  when  the  mean  strain  is  0.3,  and  less  than  6.6%  when  the  strain 
is  0.7  .  We  observe  the  formation  of  the  shear  band. 


Dynamic  case: 

•  Contrary  to  the  quasi-static  deformation,  the  beginning  of  the  localization  in  the 
dynamic  case  occurs  at  lower  shear  strains.  A  distinct  shear  band  appears  in  the 
middle  part  of  the  shear  zone  with  inclination  <|>  =  51°  to  the  x1  axis.  In  addition, 
an  intense  shear  strain  is  observed  in  the  specimen  and  the  support  areas  of  the 
shear  device. 

•  Numerical  calculations  were  not  conducted  for  advance  localization,  when  the 
results  are  mesh-sensitive. 


•  The  performed  numerical  calculations  enabled  us  to  evaluate  the  optimal 
dimensions  of  the  double  shear  specimen. 


Industrial  Fracture  Consortium  Meeting,  MIT  Cambridge,  October  13-14,  2005 


35 


Institute  of  Fundamental  Technological  Research 


Tension  device  in  testing  machine 


A  new  tension  device  was  used  by  us  to  perform 
tests  at  high  strain  rates  on  specimens  made  of 
metal  sheets. 


The  loading  and  the  displacements  of  this  device  are 
controlled  by  SHPB  acting  in  compression. 


The  role  of  the  special  device  is  to  transform  the 
compression  into  tension  (inversion  of  motion). 


The  mechanical  impedance  of  two  parts  of  the  shear 
device  and  SHPB  are  the  same. 
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Comparison  of  simple  shear  test  with  tension  test 


Equivalent  stress  [MPa] 


The  comparison  of  the  result  obtained  in  the  simple 
shear  test  and  tension  test  for  ES  steel  at  strain  rate 
100  s'1.  We  assume  the  Huber-Mises-Hencky  yield 
criterion:  y  =  yf3sn  and  <Jn=^crn  . 


The  influence  of  stress  components  &11  and  a22  must  be  taking  account. 
The  stress  intensity  is  defined  as: 

<J i  —  [<JU  +  CT 22  ^ \  1  ^22  ] 
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CONCLUSIONS 

•Considerable  homogeneity  of  the  permanent  strain  field  at  finite 
deformations  over  the  total  length  of  the  specimens  is  observed  in 
experiments  of  simple  shear  and  in  the  results  of  simulation. 

•The  proposed  method  is  the  only  known  test  providing,  in  the  case  of  a 
thin  sheet,  homogeneous  stress  and  strain  fields  in  both  the  dynamic  and 
static  tests. 

•Simple  shear  test  is  particularly  attractive,  since  the  application  of  this 
type  of  loading  path  can  result  in  large  strains  without  the  occurrence  of 
plastic  instability. 
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